Search for: All records

A new version of the US National Science Foundation National Center forAtmospheric Research (NSF NCAR) thermosphere-ionosphere-electrodynamicsgeneral circulation model (TIEGCM) has been developed and released. Thispaper describes the changes and improvements of the new version 3.0since its last major release (2.0) in 2016. These include: 1) increasingthe model resolution in both the horizontal and vertical dimensions, aswell as the ionospheric dynamo solver; 2) upward extension of the modelupper boundary to enable more accurate simulations of the topsideionosphere and neutral density in the lower exosphere; 3) improvedparameterization for thermal electron heating rate; 4) resolvingtransport of minor species N(2D); 5) treating helium as a major species;6) parameterization for additional physical processes, such as SAPS andelectrojet turbulent heating; 7) including parallel ion drag in theneutral momentum equation; 8) nudging of prognostic fields near thelower boundary from external data; 9) modification to the NO reactionrate and auroral heating rate; 10) outputs of diagnostic analysis termsof the equations; 11) new functionalities enabling model simulations ofcertain recurrent phenomena, such as solar flares and eclipse. Wepresent examples of the model validation during a moderate storm andcompare simulation results by turning on/off new functionalities todemonstrate the related new model capabilities. Furthermore, the modelis upgraded to comply with the new computer software environment at NSFNCAR for easy installation and run setup and with new visualizationtools. Finally, the model limitations and future development plans arediscussed. 

Abstract Motivated by numerous lower atmosphere climate model hindcast simulations, we performed simulations of the Earth's atmosphere from the surface up through the thermosphere‐ionosphere to reveal for the first time the century scale changes in the upper atmosphere from the 1920s through the 2010s using the Whole Atmosphere Community Climate Model—eXtended (WACCM‐X v. 2.1). We impose solar minimum conditions to get a clear indication of the effects of the long‐term forcing from greenhouse gas increases and changes of the Earth's magnetic field and to avoid the requirement for careful removal of the 11‐year solar cycle as in some previous studies using observations and models. These previous studies have shown greenhouse gas effects in the upper atmosphere but what has been missing is the time evolution with actual greenhouse gas increases throughout the last century, including the period of less than 5% increase prior to the space age and the transition to the over 25% increase in the latter half of the 20th century. Neutral temperature, density, and ionosphere changes are close to those reported in previous studies. Also, we find high correlation between the continuous carbon dioxide rate of change over this past century and that of temperature in the thermosphere and the ionosphere, attributed to the shorter adjustment time of the upper atmosphere to greenhouse gas changes relative to the longer time in the lower atmosphere. Consequently, WACCM‐X future scenario projections can provide valuable insight in the entire atmosphere of future greenhouse gas effects and mitigation efforts. 

Abstract A companion paper by Lund et al. (2020) employed a compressible model to describe the evolution of mountain waves arising due to increasing flow with time over the Southern Andes, their breaking, secondary gravity waves and acoustic waves arising from these dynamics, and their local responses. This paper describes the mountain wave, secondary gravity wave, and acoustic wave vertical fluxes of horizontal momentum, and the local and large-scale three-dimensional responses to gravity breaking and wave/mean-flow interactions accompanying this event. Mountain wave and secondary gravity wave momentum fluxes and deposition vary strongly in space and time due to variable large-scale winds and spatially-localized mountain wave and secondary gravity wave responses. Mountain wave instabilities accompanying breaking induce strong, local, largely-zonal forcing. Secondary gravity waves arising from mountain wave breaking also interact strongly with large-scale winds at altitudes above ~80km. Together, these mountain wave and secondary gravity wave interactions reveal systematic gravity-wave/mean-flow interactions having implications for both mean and tidal forcing and feedbacks. Acoustic waves likewise achieve large momentum fluxes, but typically imply significant responses only at much higher altitudes. 

Abstract This paper addresses the compressible nonlinear dynamics accompanying increasing mountain wave (MW) forcing over the southern Andes and propagation into the mesosphere and lower thermosphere (MLT) under winter conditions. A stretched grid provides very high resolution of the MW dynamics in a large computational domain. A slow increase of cross-mountain winds enables MWs to initially break in the mesosphere and extend to lower and higher altitudes thereafter. MW structure and breaking is strongly modulated by static mean and semidiurnal tide fields exhibiting a critical level at ~114 km for zonal MW propagation. Varying vertical group velocities for different zonal wavelengths λ x yield initial breaking in the lee of the major Andes peaks for λ x ~ 50 km, and extending significantly upstream for larger λ x approaching the critical level at later times. The localized extent of the Andes terrain in latitude leads to “ship wave” responses above the individual peaks at earlier times, and a much larger ship-wave response at 100 km and above as the larger-scale MWs achieve large amplitudes. Other responses above regions of MW breaking include large-scale secondary gravity waves and acoustic waves that achieve very large amplitudes extending well into the thermosphere. MW breaking also causes momentum deposition that yields local decelerations initially, which merge and extend horizontally thereafter and persist throughout the event. Companion papers examine the associated momentum fluxes, mean-flow evolution, gravity wave–tidal interactions, and the MW instability dynamics and sources of secondary gravity waves and acoustic waves. 

The mesospheric polar vortex (MPV) plays a critical role in coupling the atmosphere-ionosphere system, so its accurate simulation is imperative for robust predictions of the thermosphere and ionosphere. While the stratospheric polar vortex is widely understood and characterized, the mesospheric polar vortex is much less well-known and observed, a short-coming that must be addressed to improve predictability of the ionosphere. The winter MPV facilitates top-down coupling via the communication of high energy particle precipitation effects from the thermosphere down to the stratosphere, though the details of this mechanism are poorly understood. Coupling from the bottom-up involves gravity waves (GWs), planetary waves (PWs), and tidal interactions that are distinctly different and important during weak vs. strong vortex states, and yet remain poorly understood as well. Moreover, generation and modulation of GWs by the large wind shears at the vortex edge contribute to the generation of traveling atmospheric disturbances and traveling ionospheric disturbances. Unfortunately, representation of the MPV is generally not accurate in state-of-the-art general circulation models, even when compared to the limited observational data available. Models substantially underestimate eastward momentum at the top of the MPV, which limits the ability to predict upward effects in the thermosphere. The zonal wind bias responsible for this missing momentum in models has been attributed to deficiencies in the treatment of GWs and to an inaccurate representation of the high-latitude dynamics. In the coming decade, simulations of the MPV must be improved. 

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to performin-situmeasurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either throughin-situmeasurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR’s Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of anin situmission in regions of interest and for various conditions of solar and geomagnetic activity. 6) Calculations of the statistical significance of obtaining the primary and derived products throughout anin situmission’s lifetime. 7) Routines for the visualization of 3D datasets of GCMs and of measurements along orbit. Daedalus MASE code is accompanied by a set of Jupyter Notebooks, incorporating all required theory, references, codes and plotting in a user-friendly environment. Daedalus MASE is developed and maintained at the Department for Electrical and Computer Engineering of the Democritus University of Thrace, with key contributions from several partner institutions. 

Abstract Fritts, Wang, Lund, and Thorpe (2022,https://doi.org/10.1017/jfm.2021.1085) and Fritts, Wang, Thorpe, and Lund (2022,https://doi.org/10.1017/jfm.2021.1086) described a 3‐dimensional direct numerical simulation of interacting Kelvin‐Helmholtz instability (KHI) billows and resulting tube and knot (T&K) dynamics that arise at a stratified shear layer defined by an idealized, large‐amplitude inertia‐gravity wave. Using similar initial conditions, we performed a high‐resolution compressible simulation to explore the emission of GWs by these dynamics. The simulation confirms that such shear can induce strong KHI with large horizontal scales and billow depths that readily emit GWs having high frequencies, small horizontal wavelengths, and large vertical group velocities. The density‐weighted amplitudes of GWs reveal “fishbone” structures in vertical cross sections above and below the KHI source. Our results reveal that KHI, and their associated T&K dynamics, may be an important additional source of high‐frequency, small‐scale GWs at higher altitudes. 

Abstract Gravity waves (GWs) and their associated multi‐scale dynamics are known to play fundamental roles in energy and momentum transport and deposition processes throughout the atmosphere. We describe an initial machine learning model—the Compressible Atmosphere Model Network (CAM‐Net). CAM‐Net is trained on high‐resolution simulations by the state‐of‐the‐art model Complex Geometry Compressible Atmosphere Model (CGCAM). Two initial applications to a Kelvin‐Helmholtz instability source and mountain wave generation, propagation, breaking, and Secondary GW (SGW) generation in two wind environments are described here. Results show that CAM‐Net can capture the key 2‐D dynamics modeled by CGCAM with high precision. Spectral characteristics of primary and SGWs estimated by CAM‐Net agree well with those from CGCAM. Our results show that CAM‐Net can achieve a several order‐of‐magnitude acceleration relative to CGCAM without sacrificing accuracy and suggests a potential for machine learning to enable efficient and accurate descriptions of primary and secondary GWs in global atmospheric models. 

Abstract. The lower-thermosphere–ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.